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Abstract 

We discuss conceptual aspects of neutrino oscillations with the main em- 
phasis on the field-theoretical approach. This approach includes the neutrino 
source and detector processes and allows to obtain the neutrino transition or 
survival probabilities as cross sections derived from the Feynman diagram of 
the combined source - detection process. In this context, the neutrinos which 
are supposed to oscillate appear as propagators of the neutrino mass eigen- 
fields, connecting the source and detection processes. We consider also the 
question why the canonical neutrino oscillation formula is so robust against 
corrections and discuss the nature of the oscillating neutrino state emerging 
in the field-theoretical approach. 



*Talk given by W. Grimus at Neutrino Mixing, Meeting in Honour of Samoil Bilenky's 70th 
Birthday, Torino, March 25-27, 1999. 
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I. INTRODUCTION 



A. History 



Neutrino oscillations play a central role in neutrino physics. The idea of neutrino os- 
cillations in analogy with K^K^ oscillations was first discussed by B. Pontecorvo in the 
late fifties The concept of 2-neutrino flavour mixing (1962) and the confirmation of 
the existence of different neutrino flavours Ve and in the same year opened the way 
to neutrino flavour oscillations 0. For a general number of neutrinos, neutrino mixing is 
denoted as 



Via = ^ Uaji^Lj with a = 6, /i, T, . . . , j = 1, 2, 3, 



(1) 



labelling neutrinos flavours (types) and mass eigenfields, respectively. 

The first theory of 2-neutrino oscillations was developed by Gribov and Pontecorvo [Q. 
The 2-neutrino oscillation probabilities in the form as they are used nowadays were formu- 
lated by Bilenky and Pontecorvo in Ref. ^ and the extension to an arbitrary number of 
neutrinos (with U real) can be found in Refs. The most general form of transition 

(a 7^ (3) and survival (a = (3) probabilities (see, e.g., Refs. P,CT) reads 



exp 



i — — 



(2) 



where U denotes the unitary mixing matrix, L the distance between source and detector and 
El, the neutrino energy. The neutrino masses rrij are associated with the mass eigenfields 
Uj. Note that the above formula assumes ultra- relativistic neutrinos. With the ordering 



nil < ^2 < 



Eq. 



is rewritten as 



'^uc-^UjiiL/ Ey 



^ |f/ft-r|f/.,f + 2 J] Re [u;^U^,Up,U:,e^Y^ (-z2vrL/L°-)] 



(3) 



where the oscillation lengths are defined by 



Am2. 



(4) 



It was a very important discovery that the oscillation pattern as described by Eqs. (|^) and 
can get significantly modified by the effect of coherent forward scattering in background 

or its equivalent form (R|) has 



matter [|TT|. Otherwise, the canonical oscillation formula 



^For a detailed account of the history of neutrino oscillations see the contribution of S.M. Bilenky 
to this meeting ( hep-ph/9908335|) . 
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resisted - at least in terrestrial experiments - all attempts to find non-negligible corrections. 



The strong evidence for neutrino oscillations found in atmospheric neutrino experiments |T2 
(for reviews on neutrino oscillations see, e.g., Ref. |T3|) constitutes a major achievement and 
progress in neutrino physics. If in addition to the solar neutrino evidence for neutrino 
oscillations also the LSND evidence will be confirmed by future experiments, then neutrino 
physics is enriched by a sterile neutrino (or sterile neutrinos), a notion coined by Pontecorvo 
in Ref. [Q. In this case, in the summation in Eq. (^ there will be sterile neutrino types in 
addition to the flavours. For a review of the physics of sterile neutrinos see Ref. [|14|. 



B. Conceptual questions 

It has been indicated in several publications that the standard derivation of Eq. @ as 
found, e.g., in Ref. ^ raises a number of conceptual questions (see, e.g., Ref. ||T5| for a clear 
exposition). Let us consider some of these questions for a moment: 

1. Under which conditions is Eq. (^ valid? 

2. Do the neutrino mass eigenstates |z/j) have definite momenta or definite energies? Or 
are both smeared out? Note that in the usual derivation it is assumed that all \i>j) 

have the same momenta p and thus energies Ej = \Jrn?j + p^- One considers their time 

evolution given by ex]i{—iEjt) and finally makes the replacement t ^ L/c {c is the 
velocity of light). The last step is necessary to make contact with realistic experiments 
where the distance L between source and detector is known, but no measurement of 



the time t can be made 16 



3. Why is it allowed to make the replacement t L/c when wave functions with definite 
momentum have infinite extent in space? 

4. Why not making the replacement t cL/vj where Vj = p/Ej is the velocity of the 
neutrino with mass rrij? 

5. Evidently, the Fock space of neutrinos with definite mass exists. Is this also the case 
for flavour neutrino states as sometimes assumed in the literature? 

One aspect of the limitation of Eq. (^ is treated in terms of a coherence length 
|l^,20]. Some of the above questions are answered in the framework of the wave packet 



approach [0,0 (see also the review |2T| where a list of references can be found), where 
the neutrino momentum is smeared out, however, the size and form of wave packet is not 
determined in this approach and remains a subject to reasonable estimates. Let us shortly 
mention the answers to the above questions offered by the wave packet approach: 

Ad 1: One necessary condition for validity of the canonical formula (0) is formulated using 
the notion of a coherence length. Any coherence length in neutrino oscillations [|l^ is 
related to the oscillation length (|^) by ^2 



E 

rcoh ^ rose /r"\ 
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where Ei, is the mean neutrino energy and AE,^ is the energy spread of the neutrino 
beam, li L ^ L^^^ holds, the term with exp(— i27rL/L^°^) in Eq. is suppressed. 

Ad 2: Obviously, smearing out the neutrino momentum also smears out the neutrino ener- 
gies. 



Ad 3: The wave packet approach deals with the time variable either by putting ||18| t = L/c 
or by averaging [1^ over t. 

Ad 4: Such a replacement in the standard derivation of Eq. is wrong because it leads to 
Ejtj = EjL/p = pL + m'jL/p and consequently misses the factor 1/2 in the exponents 

MM- 



Ad 5: In Ref. El] arguments are given that a Fock space of flavour neutrinos does not exist. 



C. The field-theoretical approach 

The idea has been put forward to include the neutrino production and detection pro- 
cesses into the consideration of neutrino oscillations |]T3|. Such an approach can be realized 



with quantum mechanics - in which case the neutrinos with definite mass are unobserved 



intermediate states between the source and detection processes ||T5[ - or with quantum field 
theory where the massive neutrinos are represented by inner lines in a big Feynman diagram 
depicting the combined source - detection process |]25| , ^ . In the following we will discuss 
the field-theoretical treatment. The aims and hopes of such an approach are the following: 

A. The elimination of the arbitrariness associated with the wave packet approach, 

B. the description of neutrino oscillations by means of the particles in neutrino production 
and detection which are really manipulated in an experiment, 

C. a more complete and realistic description in order to find possible limitation of formula 

in specific experimental situations. 

Note that in this approach the question of a Fock space for flavour states has no relevance. 

Considering laboratory experiments, there are two typical situations for neutrino os- 
cillation experiments. The first one is decay at rest (DAR) of the neutrino source. Its 
corresponding Feynman diagram is depicted in Fig. |l]. The wave functions of the source and 
detector particles are localized (peaked) at xs and xd, respectively. The other situation is 
decay in flight (DIE) of the neutrino source as represented in Fig. ^ where it is assumed 
that a proton hits a target localized at xt- The detector particle sits again at xd but the 
source is not localized. In both situations the distance between source and detection is 
given by L = \xe> — xs\- Note that in the Feynman diagrams of Figs. and ^ the neutrinos 
with definite mass occur as inner lines. In the spirit of our approach, neutrino oscillation 
probabilities are proportional to the cross sections derived from the amplitudes represented 
by these diagrams. 

In the following we will work out the amplitude for neutrino oscillations as described 
by Figs. 1^ and |^ in the limit of macroscopic L, discuss the consequences and compare this 



4 



approach with the wave packet approach. Note that within the field-theoretical approach 
the notion coherent refers to a summation in the amplitude, whereas incoherent means a 
summation over squares of amplitudes, i.e., a summation in the cross section [p6| . 



II. ASSUMPTIONS AND THE RESULTING AMPLITUDE 
A. The basic assumptions 

The further discussion is based on the following assumptions: 

I. The wave function 0d of the detector particle does not spread with time, which 
amounts to 

0(f , t) =^d{x- xd) e"'^^^* , (6) 

where Edp is the sharp energy of the detector particle and ipDiv) is peaked aX y = 0. 

II. The detector is sensitive to momenta (energies) and possibly to observables commuting 
with momenta (charges, spin). 

III. The usual prescription for the calculation of the cross section is valid. 

With the amplitude symbolized by Fig. |I| the oscillation probability is obtained by 



i^a^ ^fi I -p J J-p ^r^Dl ^i^Dni 



A 



{-) (-) 



(7) 



In this equation we have indicated the average over some region V in the phase space of the 
final particle of the detection process. If no final particle of the neutrino production process 
is measured then one has to integrate over the total phase space of these final states. For 
DIF symbolized by Fig. ^ there is an additional integration J dPx over the final particles in 

the target process. By definition, at the source (detector) a neutrino 'i^i ('z/^) is produced 
(detected) if there is a charged lepton (/5^) among the final states. 

In perturbation theory with respect to weak interactions, according to the Feynman dia- 
grams Figs. and|^ one has to perform integrations J d'^Xi, J d^X2 and j d'^q corresponding 
to the Hamiltonian densities for neutrino production, detection and the propagators of the 
mass eigenfields. These integrations are non-trivial because ipD and the source (target) wave 
functions are not plane waves, but are localized at xd and xs {xt), respectively. 



B. Integration over d'^X2 

Let us discuss now in detail the integration J d'^X2 at the interaction vertex of the detec- 
tor process. For definiteness we consider transitions, which is the typical situation 
in experiments with DAR of the neutrino source (see LSND, KARMEN and reactor exper- 
iments). A remark is at order concerning the integration over dx\. We will perform this 
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integration from — oo to +00, which corresponds to initial time ti = —00 and final time 
tf = +00 of the combined source - detection process. This is certainly only approximately 
correct for DAR (see Fig. |^), because in this case the time evolution of the source particle 
should rather be considered starting at a finite time, after the preparation of the state of the 
source particle. We will later comment on this problem in detail. Looking at the neutrino 
line in Figs. |l| and ^ which connects the neutrino source with the detector particle, with As- 
sumption II saying that final states are described by plane waves, we read off the following 
integral: 



" y .-iq-{xi-X2) 



(2vr) 



e 



HI ■ 



3 



g2 _ ^2 _|_ Pj 



^-iq-{xD-Xi) 



(27r)3/2 



E 2 ' ' 2I • u;, + p'd) (8) 

q-^ — m, + le 
j ■' 



with the definitions 



Pd = PDb and Ed = Y^ E'^^ - Edp ■ (9) 
6=1 6=1 

The Fourier transform of ijj is denoted by ip. Note that xi is identical with the source 
position xs for DIF. One can show for DAR with a wave function of the neutrino source 
analogous to the detector wave function @ that integration over xi causes the replacement 
of Xi by Xs in the expression (W) p5|. 



C. The integration over d'^q and the asymptotic hmit L ^ 00 

The next step in treating the expression (H) is obviously the integration over d^q. Ac- 
tually, we need the result of this integration only in the asymptotic limit of macroscopic 
distance between source and detector L = \xd — xs\- This limit is handled by the following 
theorem [^]. 

Theorem: Suppose we have an integral of the form 

where $ is three times continuously differentiable and $, Vj$, Vj Vfc$ all decrease at least 
like l/(f^ for L = |L| — >■ 00. Furthermore, A is a constant, i.e., independent of q and L. 
Then in the asymptotic limit L ^ 00 the integral ([TD|) is given by 
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J(L) = -^^$|-VIy )e*^VO(L-3/2) for A>0, (11) 




L 

J{L) = {L-'') for A < . (12) 

Applying this theorem to the integration j d^qe~^'^'^ ■ ■ ■ in the expression (H), for each j 
we identify A with (g°)^ — m| and L with xd — xs- Apart from a factor, the result of the 
integration and the asymptotic limit L — oo, which picks out the term proportional to L~^, 
is described by the following operations on the right-hand side of Eq. (H) for each j: 

1. The denominator in the neutrino propagator is removed. 

2. In the three places in expression (§) where q occurs, i.e., in the exponential, in the 
numerator of the neutrino propagator and in ^/'/j, the replacement 



q —Qj^ with qj = \/ Ej^ — and i = {xn — Xs)/L (13) 



has to be made. 

This amounts to having neutrino j on mass shell with four-momentum 



Note that with kj the numerator of the neutrino propagator is trivially rewritten as 

— + TTij = — 'i^jikj, s)v{kj, s) , (15) 

±s 

which is the projection operator on the negative-energy solutions of the free Dirac equation. 
Of course, due to the left-handedness of the weak charged currents the term with negative 
helicity is negligible in Eq. ([T5|) . 

This consideration conforms to our physical intuition: if we have a production process 
for a particle with which one performs a scattering experiment at a distance macroscopi- 
cally separated in space from the production process, then the amplitude of the combined 
production - scattering (detection) process is the product of the production and scattering 
amplitudes. This is the physical content of the above theorem. In the case of neutrino 
oscillations the particles with which the scattering is performed are the massive neutrinos. 

In the case of neutrino mixing (|ID several different massive neutrinos are produced and 
the amplitude is the sum over the products of production and detection of the neutrino 
mass eigenstates with mass rrij. In this way interference effects arise if certain conditions, 
which will be explored in the following, are fulfilled. In summary, after having performed 
the integrations over d'^X2 and d'^q, in the asymptotic limit L ^ oo only the neutrinos on 
mass shell contribute to the amplitude in Eq. (0) which can be written as [p5| , p^ 



^Z^u, = E A^^fUp,U:^e^'^^ , (16) 
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where ^4.^ and Af denote the amphtudes for production and detection of a neutrino with 
mass rrij. Note that is the energy on the neutrino hne in Figs. |l| and |^. It is independent 
of rrij and determined by the energies of the final states of the detection process. This 
is an immediate consequence of the assumptions in the previous section, in particular, of 
Assumption I. Furthermore, due to the above-mentioned integrations and the asymptotic 
limit we obtain 

Af ^i^Di-qZ+PD)- (17) 



In the case of DAR an analogous function ijjs appears {ijjs is localized at xs), such that |25,26 

Aj (x-^siqj^+p's) with p's = ^P'sb- (18) 

b=l 

The amplitude (0) looks very similar to the amplitude appearing in the standard oscilla- 
tion formula (H), and the question when the standard formula arises amounts to a discussion 
of the conditions under which the amplitudes Aj and Af are independent of j. This will 
be done in the next section. 

In Ref. it has been shown that the finite lifetime of the neutrino source particle can 
be incorporated with the help of the Weisskopf-Wigner approximation, thus performing per- 
turbation theory starting at a final initial time tj = 0. Consequently, in the d'^X2 integration, 
the time variable X2 is integrated from to cx) and, therefore, the dq^ integrand is not simply 
given by 6{q^ + Ed) (see Eq. (||)). Nevertheless, it has been demonstrated in Ref. that 



in the limit of macroscopic L one recovers this delta function and all our manipulations at 
the detector vertex and the inner neutrino line discussed in this section are valid also for 
explicit finite lifetime of the neutrino source. However, there are effects stemming from the 
source vertex. These will be discussed in the next section without derivations. For details 
we refer the reader to Refs. PBip? . 



III. RESULTS 



The preceding discussion based on the assumptions stated in Section |II A| leads us to the 
following conclusions: 

(i) Since in the asymptotic limit L —>■ 00 the neutrinos on the inner line in the Feynman 
diagrams Fig. |I] and |^ are on mass shell we are allowed to make the interpretation that 
we have neutrino mass eigenstates |z/j) characterized by the energy E^, = E^ Vj and 

momenta qj = \ E'j-, — m?-. Thus they have all the same energy determined by the 



detection process, but the momenta are different [^. The same conclusion has also 



been reached in |T^J2^, however, for other reasons. 

The summation over Eo is incoherent, i.e., it occurs in the cross section (see Eq. (^), 
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not in the amplitude ([l6D. In this sense there are no neutrino wave packets^ in exper- 
iments conforming to our assumptions |26 . 



(iii) It has been shown in Ref. that the relation (|^) holds irrespective of the coherent 
or incoherent nature of the energy spread. However, with the two points just made 
it follows that any coherence length originates in an incoherent energy spread and 
reflects the inability to measure the energies of the final states of the detector process 



more precisely than AE,^ Note that in ||22| it has already been pointed out that 
a coherent or incoherent neutrino energy spread cannot be distinguished in neutrino 
oscillation experiments. Since the neutrino energy can in principle be determined 
with arbitrary precision the coherence length can theoretically be increased solely by 
detector manipulations p2| , |29[ . 



(iv) From Eq. (|T^ it follows that with Am\ - = ml — {k > j) the condition 

Attt,^ 1 
q3-qk^^<^D or a,D<^Lt^ (DAR + DIF) (19) 

is necessary for having neutrino oscillations with Am\- (see Eq. (^), where o"d and 
a^D are the widths of the wave function of the detector particle in momentum and 
coordinate space, respectively. Otherwise, the standard formula @) has to be corrected 
by suppressing the oscillatory term containing AmL-. The second part of equation 



( p!9D is obtained by rewriting the first part using the oscillation length (^) [Il8| , p!5| , |25 
In realistic experiments condition (0) holds because axu is a microscopic quantity, 
whereas, dropping from now on the indices j, fc, the oscillation length 

~ 2.48 m ( ( (20) 

is a macroscopic quantity. 



(v) In the case of DAR a condition analogous to (|T9|) 



axs<j^L°'' (DAR) (21) 

exists for the width of the neutrino source wave function. Note that the conditions 
( p!9D and (|21|) are independent of L, i.e., Eqs. ([T9|) and ( pi]) do not define a coherence 
length. Such conditions were first discussed in the wave packet approach by Kayser 



(vi) Now we want to consider effects of the finite lifetime ts = l/F of the neutrino source. 
The decay width F introduces an energy spread AE^, ~ F. However, the corresponding 



^The same conclusion but for a different reason was presented in Ref. p^ . 
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coherence length is much too long to be of interest in realistic experiments if we 
take the decay widths of muons or charged pions. If the neutrino source particle can be 
conceived as a free, unbound particle its wave functions spreads during its lifetime. It 



has been shown that, in order not to wash out neutrino oscillations, the condition 



msE_ 



D 



< r (22) 



has to be fulfilled, where as is the spread of ips- Defining a velocity spread of the 
source particle at rest by as/ms, where ms is its mass, we can rewrite Eq. (|22|) as 

AvsTs<^L^'\ (23) 

We can interpret this relation by saying that the spreading of the source wave function 
during the lifetime of the source particle should be less than the oscillation length in 
order not to destroy the oscillation pattern. 

Some comments concerning these results are at order. The points (i) - (iv) refer to only 
to the detector. In particular, (i) - (iii) mean that the nature of the neutrino energy spread 
is determined by the properties of the detector. This is an immediate consequence of the 
Assumptions I, II, III in Section [II A| . Making a simple numerical estimate of Eq. (^) for the 
LSND ig and KARMEN experiments, where the source particles are stopped muons. 



we have as < 0.01 MeV Am^ ~ 1 eV^ and ~ 30 MeV. Thus we get Am'^as/m^E^ < 
3 X 10~^^ MeV -C ^ 3 X 10~^^ MeV. Therefore, the coherence condition (^) is very 
well fulfilled (see Ref. for a detailed derivation and discussion of this condition). The 



quantum-mechanical conditions (|T9|) and (^Tj) have classical analogues: the uncertainty in 
the position of the detector (source) particle has to be less than the oscillation length, 
otherwise oscillations are washed out by these macroscopic, classical averaging processes . 
Consequently, if these classical conditions, which involve macroscopic lengths, are fulfilled, 
then the conditions (|I9|) and (|2T| ) involving the microscopic spread of wave functions hold 
to a much better degree. In summary, in the framework discussed here all corrections to 
Eq. (PP are negligible. 



IV. WHAT REMAINS TO BE DISCUSSED 

It remains to point out what is still missing or what was not discussed in the field- 
theoretical approach. In the case of DIP, conditions concerning the target process (replacing 
conditions concerning the source process for DAR) have to be worked out. For an attempt 
in this direction see Ref. |Q. Note that we have not been able to build in the interaction of 



the final state particles in the source with the environment. This is a potential source of a 
finite coherence length |T^,^ according to the following consideration: neutrino emission is 



interrupted at the time when one of the final particles in the source process interacts with 
the environment. This leads to a finite length of the neutrino wave train and thus to 
an energy spread AEi, oc l/iu, resulting in a finite coherence length (|^). Such a coherence 
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length might play a role for solar neutrinos with Am^ ~ 10~^ eV^ (MSW effect [0]), in 
which case the mass- squared difference and L = 150 x 10^ km could be large enough that the 
neutrino wave packets corresponding to the mass eigenstates are separated when they arrive 
on earth [r^j2^ . However, the different averaging processes involved in the calculation of 
the solar neutrino flux lead in any case effectively to the consideration of the neutrino state 
arriving on earth as an incoherent mixture of neutrino mass eigenstates (see, e.g., Ref. [0). 
A further problem is the inclusion of the MSW effect in the field-theoretical approach. This 
problem was studied in Ref. 



35 



In the theorem discussing the asymptotic limit L — >■ cxd 
(see Section [11 C\j we have assumed that L pertaining to a neutrino oscillation experiment is 
so large that in the specific situation L is in the asymptotic region of the integral (^. It was 
found in Ref. that under certain conditions there is an intermediate range for L which 
results in an oscillation formula different from (^. There is also the possibility that in some 
cases (e.g., the KARMEN experiment) it is not realistic to use the conventional procedure 
to calculate the cross section (0) by taking the asymptotic limit of the final time to infinity. 



V. CONCLUDING REMARKS 

The field-theoretical treatment of neutrino oscillations includes the neutrino production 
and detection processes and allows thus a description in terms of quantities which are really 
manipulated in an experiment. This treatment has demonstrated the robustness of the 
canonical oscillation formula (0) against corrections in realistic experimental situations. The 
reason is that in the field-theoretical approach, for the validity of Eq. (H), only the conditions 
(PD|), ( prp and (P^D, which do not involve L, have to be fulfilled. The first two conditions 
simply state that the macroscopic oscillation length (pOD has to be much larger than the 
microscopic spatial extension of the wave functions of the neutrino detection and source 
particles. Condition (^2]) refers to unstable neutrino sources and seems to be well fulfilled in 
view of ||30| Am^ < 2 eV^ for all realistic sources. Conditions involving L (oscillation lengths) 
stemming from interactions of final state particles of the source process with surrounding 
matter have not yet been included in the field-theoretical approach. However, we do not 
expect that their inclusion drastically changes the heuristic estimates of coherence lengths. 

A theoretical aspect of the approach presented here is that the properties of the oscillating 
neutrino state is largely determined by the detector. With the reasonable assumptions 
introduced in Section |II A| we arrive at the conclusion that no neutrino wave packets are 



needed for the description of neutrino oscillations. In some sense, neutrino wave packets are 
replaced by the wave functions of the neutrino detector and source (target) particles. Using 
neutrino flavour states 



,;t,a;) = ^f/:>,)e-^(^''*-*-^) with q, = J El - m] (24) 



and {i'j\i'k) = Sjk in the discussion of neutrino oscillations and performing neutrino energy 
summations in the probabilities (0) but not in the states (p^ comes close to the spirit of 
the field-theoretical approach. 
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FIGURES 




FIG. 1. Feynman diagram for decay at rest (DAR) of the neutrino source particle. The source 
(S) and detector (D) processes are symboUzed by the circles. The labels xs and xd represent the 
coordinates where the wave functions of the source and detector particles are peaked, respectively. 
We have also indicated the ns (no) momenta of the final particles originating from the source 
(detector) process and the neutrino propagator of the neutrino field with mass rrij. 




FIG. 2. Feynman diagram for decay in flight (DIF) of the neutrino source particle which is 
produced by a proton with momentum P hitting a target (T) particle localized at xt- In addition 
to the final momenta of DAR there are ut final momenta originating from the target process. 
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